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ABSTRACT 



Context. Observations of magnetic clouds (MCs) are consistent with the presence of flux ropes detected in the solar wind (SW) a few 
days after their expulsion from the Sun as coronal mass ejections (CMEs). 

Aims. Both the in situ observations of plasma velocity profiles and the increase of their size with solar distance show that MCs are 
typically expanding structures. The aim of this work is to derive the expansion properties of MCs in the inner heliosphere from 0.3 to 
1 AU. 

Methods. We analyze MCs observed by the two Helios spacecraft using in situ magnetic field and velocity measurements. We split 
the sample in two subsets: those MCs with a velocity profile that is significantly perturbed from the expected linear profile and those 
that are not. 

From the slope of the in situ measured bulk velocity along the Sun-Earth direction, we compute an expansion speed with respect to 
the cloud center for each of the analyzed MCs. 

Results. We analyze how the expansion speed depends on the MC size, the translation velocity, and the heliocentric distance, finding 
that all MCs in the subset of non-perturbed MCs expand with almost the same non-dimensional expansion rate (f). We find departures 
from this general rule for f only for perturbed MCs, and we interpret the departures as the consequence of a local and strong SW 
perturbation by SW fast streams, affecting the MC even inside its interior, in addition to the direct interaction region between the SW 
and the MC. We also compute the dependence of the mean total SW pressure on the solar distance and we confirm that the decrease 
of the total SW pressure with distance is the main origin of the observed MC expansion rate. We found that f was 0.91 ± 0.23 for 
non-perturbed MCs while f was 0.48 ± 0.79 for perturbed MCs, the larger spread in the last ones being due to the influence of the 
solar wind local environment conditions on the expansion. 
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1. Introduction 

Magnetic clouds (MCs) are magnetized plasma structures form- 
ing a particul ar subset of in terplanetary coronal mass ejections 
(ICMEs, e.g.. lBurlagalll995l) . MCs are transient structures in the 
solar wind (SW) defined by an enhanced magnetic field with re- 
spect to that found in the surrounding SW with a coherent rota- 
tion of the field of the order of about a day when these structures 
are observed at 1 AU dBurlaga et al.lll 98lb . A lower proton tem- 
perature than the expected one in the SW with the same velocity 
is an other signature of MCs tha t complement their i dentification 
(e.g jRichardson & Candll995UMarsch et alj|2009l) . 

MCs interact with their environment during their journey in 
the solar wind (SW) from the Sun to the outer heliosphere and, 
since the SW pressure (magnetic plus plasma) decreases for in- 
creasing heliocentric distance, an expansion of MCs is expected. 
In a heliospheric frame, the in situ observed bulk plasma veloc- 
ity typically decreases in magnitude from the front to the back 
inside MCs, confirming the expectation that MCs are expand- 



ing objects in the SW. Furthermore, from observations of large 
samples of MCs observed at different heliocentric distance, it has 
been sho wn that the size of M Cs increases for larger heliocentric 
distances dLeitner et al.ll2007l and references therein). 

These structures have an initial expansion from their origin 
in the Sun, as shown from observations of radial expansion at the 
corona; e.g., an example of the leading ed ge of a CME travelin g 
faster than its core is shown in Figure 6 of lTripathi et al.l(l2006l) . 
However their subsequent expansion mainly will be given by 
the envi ronmental (SW) condition s as a consequence of force 
balance dDemoulin & Dassoll2009h . 

Dynamical models have been used to describe clouds 
in e xpansion, either considering only a radial expansion 



(e.g. iFarrugia et al.ll993b Osherovich et al.ll993llFarrugia et al.l 



119971: iNakwacki et all l2008bl) . or expansion in bo t h the 
radial and axial directions (e.g., IShimazu & Vandas 20021; 



Berdichevskv etaTI l2003t iDasso et al.l l2007t iNakwack ietall 



Send offprint requests to: A.M. Gulisano 



l2008at iDemoulin & Dassoll2009l) . The main aim of these mod- 
els is to take into account the evolution of the magnetic field 
as the MC crosses the spacecraft. Another goal is to correct 
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Fig. 1. Examples of two analyzed MCs that are not significantly perturbed by a fast flow. The MC center was observed at 07-Jan- 
1975 10:39 and 04-Mar-1975 21:37 UT, for panels (a,c) and (b,d), respectively. The vertical dashed lines define the MC boundaries. 
(a,b) V x is the observed velocity component in the radial direction from the Sun, expressed in km per second. The straight line is 
the linear least square fit of the velocity in the time interval where an almost linear trend is present (where the observations are 
presented as a solid line). The linear fitting is extrapolated to the borders of the MC, which are marked with circles. (c,d) B y is the 
magnetic field component, in nT, both orthogonal to the MC axis and to the spacecraft trajectory, while the solid line represents F y , 
which is the accumulated flux of B y (Eq. |7}. The extremum of F y (proxy of the cloud center) is indicated with diamonds (A color 
version is available in the electronic version). 



the effect of mixing spatial-variation/time-evolution in the one- 
point observations to obtain a better determination of the MC 
field configuration. The expansion of several magnetic clouds 
has been analyz ed previously by fitt i ng different velocity mod- 
els to the data dFarrugiaet alj|1993t IShimazu & Vandasll2002l 
Berdichevskv et al.l 12003 TVandas et al.ll2.Q05l: lYurchvshvn et al 



20061: iDasso et al] (20071: IMandrini et alJl2007t iDemoulin et al 
20081) 



The expansion of some MCs is not always well marked with 
in situ velocity measurements. This is in particular the case for 
small MCs or those overtaken by fast streams. Slow magnetic 
clouds, with velocities lower than or of the order of 400 km/s, in 
general have small sizes, low magnetic field stre ngths, and only 
a few of them present shocks or sheaths (e.g., iTsurutani et al.l 
12004 . Fast streams overtaking magnetic clouds from behind 
can compress the magnetic field in the rear for the overtaken 
MC, for instance in some cases forming large struc tures called 
merged interaction regions (e.g. jBurlagaet al.ll2003l) . The inter- 
action between a stream and an MC can affect the internal struc- 
ture of the cloud ( e.g., as shown from numerical simulations by 
IXiong et"a"Dl2007l) . The difference between the velocities of the 
front and back boundaries, called AV b s , was frequently used to 
qualify how important the expansion of an observed MC is. A 
larger AV b s is favorable for the presence of shocks surround- 



ing the MC, especia lly for the presence of a backward shock 
dGosling et al.l Il994h . A large AV b s is less important for the 
presence of a frontal shock since a frontal shock is also created 
by a large difference between the MC global velocity and the 
overtaken SW velocity. 

The quantity AV b s is a good proxy of the time variation of 
the global size of the MC, however, AV b s does not express how 
fast the expansion of an element of fluid is, since AV b s depends 
strongly on how big the studied MC is. For e xample the MC 
obser ved by ACE at 1 AU on 29 October 2004 (IMandrini et al.l 
120071) is formed by a flux rope with a large radius, i? a; 0.17 AU, 
and it also has a large AV b s ~ 400 km s -1 , and so, at first sight, 
it can be qualified as a very rapidly expanding MC. However, let 
suppose that the same MC would have most its flux having been 
reconnected with the encountered S W during t he transit from the 
Sun, a s has been observed in some cases (e.g jDasso et al 120061 
120071) . so that only the flux rope core would have been observed 
as a MC. If the remaining flux rope would have a radius of only 
10~ 2 AU, it would have shown AV obs ~ 400/17 a 24 km s" 1 , so 
it would have been qualified as a slowly expanding MC. 

More generally, small flux ropes are expected to have in- 
trinsically small AV b s , an expectation confirmed by the data 
(Fig.|3^,b). MCs have a broad range of sizes, with flux rope radii 
of a few 0.1 AU down to a few 10 3 AU dLvnch et al.l 12003b 
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iFeng etal]l2007h . and it is necessary to quantify their expan- 
sion rate independently of their size. In this study, we analyze 
the expansion of MCs in the inner heliosphere, and find a non- 
dimensional expansion coefficient (£), which can be quantified 
from one-point in situ observations of the bulk velocity time pro- 
file of the cloud. We demonstrate that f characterizes the expan- 
sion rate of the MC, independently of its size. 

We first describe the data used, and then the method to define 
the main properties of the MC (Sect. |2). In Sect. [3] we analyze 
the properties of the MC expansion, defining a proper expansion 
coefficient. We derive specific properties of two groups of MCs, 
defined from their interaction with the SW environment. Then, 
we relate the MC expansion rate to the decrease of the total SW 
pressure with solar distance. We summarize our results in Sect.|4] 
and conclude in Sect. [5] 

2. Data and method 

2.1. Helios data base 

We have studied the MCs reported by different authors from the 
Helios 1 and 2 missions (Bothmer& Schwennl[l998l: iLiu et alJ 
120051: iLeitneret all 120071) : from November 1974 to 1985 for 
Helios 1 and from January 1976 to 1980 for Helios 2. We 
analy zed observations of plasma properties dRosenbauer et al.1 
|1977 | ). in particula r bulk velocity and dens ity of protons 
(iMarsch et all 19821) . and magnetic field vector (INeubauer et al.l 
119771) . for a time series with a temporal cadence of 40 ± 1 sec- 
onds. 

The magnetic and velocity fields observations are in a right- 
handed system of coordinates (x, y, z). x corresponds to the Sun- 
Spacecraft direction, y is on the ecliptic plane and points from 
East to West (in the same direction as the planetary motion), 
and z points to the North (perpendicular to the ecliptic plane and 
closing the right-handed system). 

2.2. Definition of the MC local frame 

To facilitate the understanding of MC properties, we define a 
system of coordinates linked to the cloud in which z c i ou d is along 
the cloud axis (with B ZjC ioud > at the MC axis). Since the ve- 
locity of an MC is nearly in the Sun-spacecraft direction and as 
its speed is much higher than the spacecraft speed (which can be 
supposed to be at rest during the cloud observing time), we as- 
sume a rectilinear spacecraft trajectory in the cloud frame. The 
trajectory defines a direction d (pointing toward the Sun); then, 
we define y c ioud in the direction z c i ou d X d and i c ioud completes 
the right-handed orthonormal base (xdoud,ycioud,Zcioud)- We also 
define the impact parameter, p, as the minimum distance from 
the spacecraft to the cloud axis. 

The observed magnetic field in an MC can be expressed in 
this local frame, transforming the observed components (B x , B y , 
B z ) with a rotation matrix to (B VjC ioud, #y,cioud> #z,cioud)- In partic- 
ular, for p = and an MC described by a cylindrical magnetic 
configuration, i.e. B(r) = B z (r)z + B^,{r)^>, we have i c i ou d = r and 
Jcioud = 4> when the spacecraft leaves the cloud. In this particular 
case, the magnetic field data will show: Becloud = 0, a large and 
coherent variation of Becloud (with a change of sign), and an in- 
termediate and coherent variation of Becloud, from low values at 
one cloud edge, with the largest value at its axis and returning to 
low values at the other edge. 

More generally, the local system of coordinates is especially 
useful when p is small compared to the MC radius (R) since the 
direction of the MC axis can be found using a fitting method or 



applying the minimum variance (MV) technique to the normal- 
ized t ime series of the observed magnetic field (e.g. lDasso et al.l 
120061 and references therein). In p articular, from the anal ysis of 
a set of cylindrical synthetic MCs, iGulisano et al.l (120071) found 
that the normalized MV technique provides a deviation of the 
real orientation of the main MC axis of less than 10° even for p 
as large as 50% of the MC radius. 

2.3. Definition of the MC boundaries 

As the first step of an iterative process, we choose the MC 
boundaries reported in the literature, and perform a normal- 
ized minimum variance analysis to find the local frame of the 
MC. We then analyze the magnetic field components in the lo- 
cal frame, and redefine the boundaries of each MC, according 
to the expected typical behavior of the axial field (B zc i ouc i, hav- 
ing its maximum near the center and decreasing toward the MC 
boundaries), the azimuthal field {B y c i ouc i, maximum at one of the 
borders, minimum at the other one and changing its signs near 
the cloud center), and B x . ^ mu j, which is expect ed to be small and 
with small variations (see lGulisano et al.l2007l) . Moreover, in the 
MC frame it is easier to differentiate the SW and MC sheath 
magnetic field, which is fluctuating, from the MC field, which 
has a coherent and expected behaviour for the three field com- 
ponents. We then moved the borders in order to reach these prop- 
erties of the local field components, and performed the same pro- 
cedure iteratively to find an improved orientation, using several 
times the minimum variance technique. We applied this proce- 
dure to each MC of our sample. 

2.4. Characterization of the MC expansion 

Most MCs have a higher velocity in their front than in their 
back, showing that they are expanding magnetic structures in 
the SW. About half of the studied MCs have well defined linear 
profile V x (t) (Fig.[TJi,b), while for the other half, V x (t) is nearly 
linear only in a part of the MC which includes the MC center 
(Fig. Hk,b). The distortions of V x (t) are more frequently due to 
an overtaking faster SW flow in the back of the MC. 

We split the data set in two groups: non-perturbed MCs for 
cases where the velocity profile presents a linear trend in more 
than 75% of the full size of the MC and perturbed MCs for cases 
where this is not satisfied. There are almost as many perturbed 
as non-perturbed MCs, considering data from each spacecraft 
separately and both of them together. The measured temporal 
profile V x (t ) is fitted using a least square fit with a linear function 
of time, 

V xM (t) = V oM + (dV x /dt) 6t t, (1) 

where (dV A -/df)fit is the fitted slope of the linear function. We 
always keep the fitting range inside the MC, but restrict it to the 
most linear part of the observed profile. This choice minimizes 
the effect of the interacting flows (but it does not fully remove it, 
since a long term interaction can, a priori, change the expansion 
rate of the full MC). 

The linear fit is used to define the velocities V x fit(Un) and 
V x fit(t out ) at the MC boundaries (Sect. 12.3b . Then, we define the 
full expansion velocity of an MC as: 

AV, = V xfit (tin) - V xM (t m „) . (2) 

For non-perturbed MCs, AV A , is very close to the observed ve- 
locity difference V A (f,„) - V x (t out ), see e.g. upper panels of Fig. [T] 
For perturbed MCs, this procedure minimizes the effects of the 
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Fig. 2. Examples of two analyzed MCs that are perturbed by a fast flow as seen in the upper panels. The MC center was observed 
at 30-Jan-1977 03:18 and 23-Jun-1980 12:25 UT, for panels (a,c) and (b,d), respectively. The external perturbations enter in a 
significant part of the MCs. The same quantities, as in Fig. Q] are shown (A color version is available in the electronic version). 



perturbations entering the MC. Then, the expansion velocity is 
defined consistently for the full set of MCs. 

2.5. MC size 

From the determination of the boundaries described above, we 
can estimate the size, S, of the flux rope along x (the Sun- 
Spacecraft direction). S is computed as the time duration of ob- 
servation of the MC multiplied by the MC velocity at the center 
of the structure (see end of section \2.7\ . 

We perform a linear least square fit in log-log plots of S as 
a function of the distance to the Sun (D). As expected, the size 
has a clear dependence with D: 



S Ml set = (0.23 + 0.01) D 

S non-perturbed = (0.32 + 0.02) D 
perturbed — 

(0.16 + 0.01) D 



0.78±0.12 



0.89±0.15 



0.45±0.16 



(3) 



where S and D are in AU. Our results are compatible within the 
error bars with previous studies: 



^Bothmer = (0.24 + 0.01) D 

S Liu = (0.25 + 0.01) D 



0.78±0.10 



,0.92±0.07 



(4) 



S Wa 



Wans 



0.19 



D 



^Leitner — 

(0.20 + 0.02) D 



,0.61 ±0.09 



where d is an est imation of the true diameter o f the MC. Our 
results are closer to iBothmer & Schwenn l dT998l and iLiu et al.l 
d2005l) who analyzed MCs and ICMEs, respectively. The larger 



difference is between our results and the last two ones while 
they a re based on the most different sets, as follows. IWang et all 
d2005h studied a large set of ICMEs defined only by a mea- 
sured temperature lower by a factor of 2 than expecte d in the SW 
with the same speed (e.g jRichardson & Cand2004l) . This set in- 
cluded MCs, but it is dominated by non-MC ICMEs. Conversely, 
iLeitner et al.l (j2007) analyzed only MCs, with a strict classical 
definition. They fitted the magnetic field observations with a 
classical cylindrical linear force-free field, then they found the 
impact parameter and the orientation of the MC axis to estimate 
the true di ameter, d, o f the MCs. So the selected events and the 
method of IWang et alj andlLeitner et al.l are noticeably different. 
We also note that lLeitner et al.1 (120071) found a larger exponent, 
1.14 + 0.44, when they restrict their data to D < 1 AU. This indi- 
cates that the relation is not strictly a power-law and this could be 
the main origin of the different results (which are so dependent 
on the distribution of events with solar distance in the selected 
sets). 



2.6. Magnetic field strength 

Another important characteristic of MCs is their magnetic field 
strength. We define the average field < B > within each MC and 
proceed as above with S . < B > has an even stronger dependence 
upon distance: 



< B > m set = (10.9 ± 0.4) D 

< B > non -perturbed — 

(11.4 ± 0.5) D 

< B > p erturbed = (10.4 ± 0.6) D 



-1.85±0.07 



1.85±0.11 



1.89±0.10 



(5) 
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Fig. 3. The panels a-d show the correlation analysis between proxies for MC expansion with different physical quantities. The MCs 
are separated in two groups: perturbed (empty square symbol) and non-perturbed (filled circle symbol). The straight lines are the 
result of a least square fit for perturbed (dashed line), non-perturbed (thin continuous line), and for both set of MCs in the list of 
events shown in Table Q] (thick continuous line). AV A is defined by Eq. ©, At = f out - t m , V c is the velocity of the MC center (or 
at the closest approach distance), and D is the distance to the Sun. The fitted values and the obtained correlation coefficient (r) 
are included as insets, considering the different groups. More significant correlation is present for the non-perturbed cases (A color 
version is available in the electronic version). 



where the units of < B > and D are nT and AU respectively. Our 
results have a stronger dependence on D than previous results: 



< B > Vl u 

< B >Wang 
< B >Leitner 



(7.4 + 0.4) D 
8.3 D 
(19+ 1.4) D 



-1.40±0.08 



-1.52 



-1.30±0.09 



(6) 



Again the s trongest diffe r ence exists between our results and 
the ones of iLeitner et al.l (120071) . The origin of this difference 
is exp ected to be the sam e as for the size. Indeed for D < 
1 AU, ILeitner et al.l d2007l) found a more negative exponent, 
- 1 .64 + 0.4, closer to our results, as ocurred for the size. 

The typical expan sion speed in MCs is o f the order of half 
the Alfven speed (e.g jKlein & Burlagall982h . In our studied set 
of MCs, we have also found that the expansion speed was lower 
than the Alfven speed [not shown]. This is a necessary condition 
to expect that the magnetic field evolves globally, adapting its 
initial magnetic field during its expansion, because the Alfven 
speed is the velocity of propagation of information through the 
magnetic structure. 



2.7. MC center and translation velocity 

Following Dasso et al. I J2006T) . we define the accumulative flux 
per unit length L (along the MC axial direction): 



F y (h,t 2 ) 



= [ By e 
Jt, 



;cloud(t) V x ,cloud(t') df . 



(7) 



Here we neglect the evolution of the magnetic field during 
the spacecraft cro ssing period (so also the "aging" effect, see 
iDasso et"aT]l2007h . The set of field lines, passing at the posi- 
tion of the spacecraft at t\, with the hyphotesis of symmetry of 
translation along the main axis, defines a magnetic flux surface, 
which is wrapped around the flux rope axis. Then, any mag- 
netic flux surface will be crossed at least twice by the space- 
craft, once at t\ and once at t 2 defined by Fy(ti,t 2 ) = 0. Then, 
this property of F y (t\ , t 2 ) permits us to associate any out-bound 
position, within the flux rope, to its in-bound position belong- 
ing to the same magnetic-flux surface. The global extremum of 
Fy(t\,t), for t\ having a fixed value, locates the position where 
the spacecraft trajectory has the closest approach distance to the 
MC axis (MC center). This position can also be found directly 
from B Y C i OU! j(t), where B >%c [ ou d crosses zero. Nevertheless, using 
the integral quantity F y has the advantage of decreasing the fluc- 
tuations and of outlining the global extremum, compared to the 
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Fig. 4. Panels a-d show the correlation analysis that tests for the dependence of the non-dimensional expansion factor £ (Eq.[8]) as a 
function of other parameters. Perturbed and non-perturbed MCs are represented as in Fig. [3] S corrected to l au and < B > C otreeted to l au 
are normalized to 1 AU using the size and field strength dependence on the distance, according to the relationship given in Eq. ((9). 
(A color version is available in the electronic version). 



local extrema (see lower panels of Figs.Q]and[2]i. 
The velocity at the MC center (V c ) is computed from the fitted 
linear regression Eq. ([1} evaluated at the time when the space- 
craft reaches the MC center. 



3. Expansion rate of MCs 

3.1. Correlation involving the expansion velocity 

From here on, we classify the MCs belonging to the full set of 
events, according to the quality of their velocity and magnetic 
observations. If they were too noisy or with a lot of data gaps, 
we exclude them from the following study, keeping only those 
MCs with relatively good quality (listed in Table [TJ. 

AV X , as defined by Eq. (O, characterizes the expansion speed 
of the crossed MC. However, as outlined in Sect. \T\ AV X is ex- 
pected to be strongly correlated with the MC size, so that it does 
not express directly how fast a given parcel of plasma is expand- 
ing in the MC. We therefore define below, after a few steps, a 
better measure of the expansion rate. The size of an MC is pro- 
portional to At = fout^fin and to V c . Fig. [3^ shows a clear positive 
correlation between AV X and At. Moreover the least square fit of 
a straight line for the full set of MCs gives a fitted curve passing 
in the vicinity of the origin (within the uncertainties present on 
the slope). This affine correlation is then removed by computing 
AV x /At. This quantity also shows, as expected, a positive corre- 
lation with V c (Fig. [3J5), but differently above, the fitted straight 
line stays far from the origin, so we cannot simply remove the 



correlation by dividing AV x /At by V c . However, its dependence 
on V% brings the fitted straight line close enough to the origin 
(within the uncertainties of the fit, Fig. [3};) so that AV x /(AtV^) 
is a meaningful quantity. These correlations are present for both 
groups of MCs, but they are much stronger for non-perturbed 
MCs (Fig. [3]5,c). Other correlations have being attempted with 
the above methodology. Either there is no significant correla- 
tion, or the fitted curve passes far from the origin. There is still 
the exception that AV x /(AtV^) has an affine correlation with 1 /D 
for non-perturbed MCs (Fig. |3jl). 

3.2. Non-dimensional expansion rate 

The above empirical correlation analysis suggests that we can 
define the non-dimensional expansion rate as the quantity: 



The first steps were to remove the MC size dependence, while 
the last step could be further justified by the need to have a non- 
dimensional coefficient. Finally, it is remarkable that the corre- 
lation analysis of the MC data leads to the d efinition of the same 
variab le, as the theoretical analysis of iDemoulin & Dassol 
J2OO9I) . 

We next verify that £ is no longer dependent on At, V c , D 
or some combination of them. Figures |4^,b show two examples 
of this exploration. Indeed, the non-perturbed MCs show almost 
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Table 1. List of MC events. 



s/c 


T c 


Group 


&VJ{t OUJ - t in 


) i 




d-m-y h:m (UT) 




i - 1 i - 1 
km s h 




HI 


07-Jan-1975 10:39 


N 


11.1 


1.0 


HI 


04-Mar-1975 21:37 


N 


8.75 


0.63 


HI 


02-Apr-1975 09:00 


P 
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The first column in Table 1 indicates the spacecraft (Helios 1 or Helios 
2), T c in the second column is the time for the observation of the MC 
center (or for the minimum approach distance). The MCs are separated 
in two groups: perturbed (P) and non-perturbed (N) by a fast S W stream. 
AV x /(t out - f,„) is minus the fitted slope of the temporal velocity profile 
(see, Eq.[2j> and ( is the non-dimensional expansion coefficient (Eq.[8j. 
kV x l{t OM -t in ) < means observed compression of the MC. The average 
values and the standard deviations are given at the bottom. 



2- 



i;=5.8 x 10 


4 AV x +0.8 ; 


r= 0.3 


(Non-perturbed) 






x 10 


3 AV x +0.092; 


r= 0.7E 


(Perturbed) 






- l,—2. 9x10 


^AV +0 34 ' 












• 








• 




B 




• 








B" " 










* □ 

□ 



-100 



100 200 

AV (km/s) 



300 



400 



500 



Fig. 5. Perturbed and non-perturbed MCs have a remarkably dif- 
ferent behavior of £ when they are plotted as a function of AV X . 
The drawing convention is the same as in Figs. I3I4K A color ver- 
sion is available in the electronic version). 




0.4 0.6 

D (AU) 

Fig. 6. Cartoon of a possible evolution of the size of the MC with 
the helio-distance, showing the expected global expansion (thick 
solid line), an example of a non-perturbed MC (thin solid) and a 
perturbed MC (dashed line). (A color version is available in the 
electronic version). 



*^ B ^corrected I 



= <B> D 



+ 1.8 



(9) 



We find that there is no significant correlation between ^ and 

S corrected to 1 AU. as Well as between ^ and < B >corrected to 1 AU 

for 

the non-perturbed MCs (Fig. @J;,d). Taking other exponents, in 
the vicinity (±0.4) of the exponents used in Eq.|9] we also find 
very low correlation coefficients, in the ranges [0.08,0.15] for 

S corrected to 1 AU and [0.15,0.29] for < B > CO rrected to 1 AU- 



no correlation, while there are still some correlations when the 
perturbed MCs are considered. Figure |4j) also shows that even 
for slow MCs (see Sect. Q]i, the non-dimensional expansion rate 
(£) is not correlated with V c . 

Still, does ( depend on the properties of the MC considered? 
To test this we first need to remove the distance dependence on 
S and < B >, found in Sect. l2.5l and l2.6l by defining values at a 
given solar distance (here taken at 1 AU). We use: 



3.3. Expansion of non-perturbed MCs 

Perturbed and non-perturbed MCs have a remarkable different 
behavior of ( as a function of AV X (Fig. [5]). While for non- 
perturbed MCs the correlations have been removed, the per- 
turbed ones have ( well correlated with AV, (r = 0.79). 

To explain the different behavior of ( let us take into account 
that the dependence of the size on the heliocentric distance is of 
the form (as observed from several statistical studies, Eqs . [3]|4j>: 



' corrected to 1 AU 



= S D 



-0.N 



S =S (D/D ) m , 



(10) 
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where So is the reference size at the distance Do. Its physical 
origin is the approxi mate pressure balance betw een the MC and 
the surrounding SW (iDemoulin & Dassoll2009l see Sect. 13.5b . 
This physical driving of the expansion is expected to induce a 
smooth expansion so that the size of an individual MC closely 
follows Eq. (TTOl i. Then, for non-perturbed MCs, we can differen- 
tiate Eq. dTOb with time to derive the expansion velocity AV X : 

dS dS dD S 

AV,~ — ~ ~ m— V c . (11) 

dt dD dt D V ; 

Then, the non-dimensional expanding rate of Eq. © is: 
_AV i D_ 

b non-perturbed — ^ ^2 ' K^^-J 

This implies that ( is independent of the size and the veloc- 
ity of the MC as well as its distance from the Sun and its global 
expansion rate AV X , in agreement with our results (Figs. 1415) . It 
also implies that the velocity measurements across a given MC 
permits to estimate the exponent m in Eq. dTOb . Indeed the val- 
ues of m deduced from the statistical study of MC size versus 
distance (Eq.O are in agreement with the mean value found in- 
dependently with the velocity measurements (Fig. [5] TableQ]). 

3.4. Expansion of perturbed MCs 

For perturbed MCs, the estimation of ( both from the size evo- 
lution with distance and from the measured velocity still gives 
consistent estimations (Eq.[3| Table[T]i. However, the mean value 
of £ is lower (£ m m ss 0.45) than for non-perturbed MCs 
(£ ^ m « 0.8). More important, ( has a much larger spread 
for perturbed than non-perturbed MCs (by a factor » 4 larger). 
Such a larger dispersion of £ is expected to be the result of the 
variable physical parameters (such as ram pressure) of the over- 
taking streams, and also because the interaction is expected to be 
in a different temporal stage for different perturbed MCs (rang- 
ing from the beginning to the end of the interaction period). 

The effect of an overtaking stream is simply to compress the 
MC (at first thought), so MCs perturbed by this effect are ex- 
pected to have a lower ( than non-perturbed MCs. This is true in 
average, but there is a significant fraction (5/16) of perturbed 
MCs that are in fact expanding faster than the mean of non- 
perturbed MCs. The largest ( is also obtained for a perturbed 
MC. Moreover, ( for perturbed MCs still has a good correlation 
with A V v , opposite to the result obtained for non-perturbed MCs 
(Fig. IS. Why do perturbed MCs have these properties? 

When an MC is overtaken by a fast SW stream, it is com- 
pressed by the ram, plasma and magnetic pressure of the overtak- 
ing stream, so its size increases less rapidly with solar distance 
(than without interaction). If the interaction is strong enough, 
this can indeed stop the natural expansion and create an MC in 
compression, as in the 3 cases present in TableQ] where AV X < 0. 
A sketch of such evolution is given in Fig. [6] However, this inter- 
action will not last for a long period of time since the overtaking 
stream can overtake the flux rope from both sides. As the total 
pressure in the back of the MC decreases, the expansion rate of 
the MC increases. Indeed, its expansion rate could be faster than 
the typical one for non-perturbed MCs, as follows. The compres- 
sion has provided an internal pressure that is stronger than the 
surrounding SW total pressure. Then, when the extra pressure 
of the overtaking stream has significantly decreased, the MC 
has an over-pressure compared to the surrounding S W, so it ex- 
pands faster than usual i. e., an overexpansion, see lGosling et al.l 
(119951) . Indeed, the flux rope is expected to evolve towards the 
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Fig. 7. Log-log plots of the S W total pressure and its components 
as a function of the solar distance. The solid line shows the least 
square straight line fitted to the total pressure as computed from 
the points marked with circles. The total pressure in the SW de- 
creases as P SW (D) = PqD 2 9 (A color version is available in the 
electronic version). 



expected size that it would have achieved without the overtaken 
SW flow. 

So, depending on which time the MC is observed in the in- 
teraction process, a perturbed MC can expand slower or faster 
than without interaction (Pig. |6j>. This explains the dispersion of 
£ found for perturbed MCs, but also the presence of some MCs 
with faster expansion than usual. 

In the case of perturbed MCs, their sizes still follow Eq. (ITQb 
on average (as shown by Eq. [3}, but it has no meaning to apply 
this law locally to a given MC. In particular, we cannot differ- 
entiate Eq. ( [TO with time to get an estimation of the local ex- 
pansion velocity of a perturbed MC (so we cannot write Eq.HTt. 
Rather we can use Eq. (TT0T > only to have an approximate size S 
in the expression of 



^perturbed — 



AV^ D_ ^ AV X D "'D 
~S~V C ~ 



1 -m 



SoV c 



(13) 



The dependence of D l ~'" « D 05 is relatively weak, but still 
present (Fig. |4^). ( also has a dependence on V c , but the range 
of V c within the studied perturbed group is very limited to de- 
rive a reliable dependence from the observations (Fig. |4j>), and, 
moreover, there is a dependence on Sq that we cannot quantify. 
It remains that the strongest correlation of ^ is with AV, (Pig. [5j. 
Indeed, for perturbed MCs, £ strongly reflects their local expan- 
sion behavior, so that ( is strongly correlated to AV X . 

3.5. Physical origin of MC expansion 

The main driver of MC expansion is t he rapid decrease of the to - 
tal SW pressure with solar distance (IDemoulin & Das"sol I2009I 
Other effects, such as the internal over-pressure, the presence of 
a shock, as well as the radial distribution and the amount of twist 
within the flux rope have a much weaker influence on the expan- 
sion. This result was obtained by solving the MHD equations 
for flux ropes having various magnetic field profiles, and with 
ideal MHD or fully relaxed states (minimizing magnetic energy 
while preserving magnetic helicity). Within the typical SW con- 
ditions, they have shown that any force-free flux rope will have 
an almost self similar expansion, so a velocity profile almost lin- 
ear with time as observed by a spacecraft crossing an MC (e.g. 
Figs. 1112b . They also relate the normalized expansion rate ( to 
the exponent tip of the total SW pressure as a function of the 
distance D to the Sun (defined by p oc D~ np ) as £ ss tip I A. 

Here we further test the above theory with the MCs analyzed 
in this paper, by comparing the value obtained for £ with the 
value of tip obtained from previous studies of the SW. 
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According to M ariam & Neubauerl (Il990l) . from fitting a 
power law to observations of the field strength in the inner helio- 
sphere according to B — Bq(D / Do)~" B , a global decay law is ob- 
tained with n B = 1.6+0.1 (B = 3.8 +0.2 nT at D = 1AU) from 
Helios 1, and n B = 1.8 + 0.1 (B = 3.3 + 0.2 nT at D = 1AU) 
from Helios 2. 

For the proton density (N p = N P $(D / Do)~" N ) we consider a 
density of iVo = 7 ± 4 cm' 3 at 1 AU (averaging slow and fast 
SW according to ISchwennll2006h and «ai = 2 (corresponding 
to the 2D expansion for the stationary SW with constant radial 
velocity). 

According to lSchwennl d2006l) and lTotten et al.l (1 19951) . it is 
possible to represent a typical dependence of the proton temper- 
ature (T p ) upon D as approximately T p = T p o(D/Do)~ nT p with 
n Tp = 1.0 + 0.1 (T pfi = (1.3 + 1.0) x 10 5 K at D = 1AU). 

For electron tem perature (T e = T e o(D/Do)~" T <) we follow 
iMarsch et al.l d 19891) . in particular their results for the ranges 
of velocities (300-500) km/s to better represent the typical 
conditions of the SW . For the velocity range (300-400) km/s, 
IMarsch et al.l ( fl989h found T efl = (1.3 + 0.4) x 10 5 K and 
n Te = 0.5 + 0.1; for the velocity range (400-500) km/s, T efi = 
(1.4 ± 0.4) x 10 5 K and n Te = 0.4 + 0.1. For electrons we then 
consider mean temperature averaged over these two ranges of 
SW speeds. 

The partial pressures (magnetic, proton, and electron) are 
shown in Fig. [7] Neglecting the small effect of the a particles, 
the total pressure in the SW (Psw) is: Psw - Pb + Pp + Pe- 
We then propose that the total pressure also follows a power law 
(Psw - PqD np ) and fit this power law to Psw(Dd- Then, we 
fit Psw with a power law (PqD^" p ). The sum of different power 
laws is generally not a power law, however in the present case we 
still find a total pressure which is very close to a power law since 
the magnetic and plasma pressures have similar exponents. This 
also implies that the exponent found, np = 2.91 +0.31, has a low 
sensitivity to the plasma B and the relative pressure contribution 
of the electrons and pr otons. 

Using the result of IDemoulin & Dassol (120091) that f « n P /4 
for force-free flux rope, we found f = 0.73 ± 0.08, in full agree- 
ment, within the error bars, with f found from velocity mea- 
surement in non-perturbed MCs (Table [1}. This further demon- 
strates that MC expansion is mainly driven by the decrease of 
the surrounding SW total pressure with solar distance. The main 
departure from this global evolution is due to the presence of 
overtaking flows. 



4. Summary and discussion of main results 

MCs have a specific magnetic configuration, forming flux ropes 
which expand in all directions, unlike the almost 2D expansion 
of the surrounding SW. But how fast do they expand? Is the ex- 
pansion rate specific for each MC or is there a common expan- 
sion rate? What is the role of the surrounding SW? Finally, what 
is the main driver of such 3D expansion? 

In order to answer these questions we have re-analyzed a sig- 
nificant set of MCs observed by both Helios spacecrafts. In order 
to better define the MC extension, we first analyzed the mag- 
netic data, finding the direction of the flux rope axis, and then 
we rotated the magnetic data in the MC local frame. This step is 
important to separate the axial and ortho-axial field components 
since they have very different spatial distributions in a flux rope, 
and since we then can use the magnetic flux conservation of the 
azimuthal comp onent as a constrain t on the boundary p ositions 
(e.g., as done in iDasso et ail d2006h : ISteed et alJ d2008l) V Then, 



in the local MC frame, we can better define the boundaries of the 
MC. 

The observed velocity profile typically has a linear varia- 
tion with time, with a larger velocity in front than in the back 
of the MC. This is a clear signature of expansion. On top of 
this linear trend, fluctuations of the velocity are relatively weak, 
with the most noticeable exception occurring when an overtak- 
ing fast stream is observed in the back of the MC. Such fast 
flow can enter the MC, removing the linear temporal trend. We 
consider these overtaken MCs in a separate group as perturbed 
MCs (Fig. [2]). We exclude from the analysis the MCs where the 
overtaking flow was extending more than half the MC size and 
MCs where the data gaps were too large. The remaining MCs 
are classified as non-perturbed, even if some of them have weak 
perturbations in their velocity profiles. These perturbations are 
filtered by considering only the major part of the velocity profile 
where the profile is almost linear with time (Fig. [1}. 

The group of non-perturbed MCs has a broad and typical 
range of sizes and magnetic field strengths (« [0.1,0.75] AU 
and x [7, 25] nT respectively when rescaled at 1 AU). They also 
have a broad range of expansion velocities (« [80, 500] km/s). 
Such a range of expansion velocity cannot be explained by the 
range of observing distances, (D in [0.3, 1] AU), since the ex- 
pansion velocity decreases only weakly within this range of D. 
However, we found that the expansion velocity is proportional 
to the MC size. By further analyzing the correlation between the 
observed expansion velocity and other measured quantities, such 
as the MC velocity, we empirically defined a non-dimensional 
expansion coefficient ( (Eq. [8]). For the non-perturbed MCs, £ 
is independent of all the other characteristics of the MCs (such 
as size and field strength). Moreover, this empirical definition 
of f , obtained by removing the correlation in the data between 
the expansion rate and other quantities, finally defines the same 
( quantity as the one d efined from theoretical considerations by 
IDemoulin et alJ (120081) . We conclude that f characterizes the ex- 
pansion rate of non-perturbed MCs. 

For the non-perturbed MCs, we found that £ is confined to 
a narrow interval: 0.91 + 0.23. This is consistent with the result 
obtained at 1 AU for a set of 26 MCs observed by Wind and ACE 
(IDemoulin et al.ll2008l) . Indeed, we found that £ is independent 
of solar distance (within [0.3, 1] AU) in the Helios MCs. 

What is the orig i n of th is common expansion rate of MCs? 
IDemoulin & Dassol d2009l) have shown theoretically that the 
main origin of MC expansion is the decrease of the total SW 
pressure with solar distance D. With a SW pressure decreas- 
ing as D~" F , they found that £ « tip/ 4 independently of the 
magnetic structure of the flux rope forming the MC. In the 
present work, we re-analyzed the total SW pressure variation 
with D, revising previous studies that also analyzed Helios data 
dMariani & Neubauer|[T990t ISchwennll2006l iTotten et alJll995t 
IMarsch et al.lll989l) . We found n P = 2.91 + 0.31, which implies 
«p/4 = 0.73 ± 0.08, in agreement with our estimation of f from 
the measured velocity in MCs. We then confirm that the fast de- 
crease of the total SW pressure with solar distance is the main 
cause of the MC expansion rate. 

For MCs overtak en by a fast SW stre am (or by another flux 
rope on its back, e.g. lDasso et al.l (120091) ). we minimize its im- 
portance in the estimation of the MC expansion rate by consid- 
ering only the part of the velocity profile which is nearly lin- 
ear with time. Still, the mean computed £ for perturbed MCs is 
significantly lower than the mean value for non-perturbed MCs, 
showing that the overtaking flows have a more global effect on 
MCs (than the part where the velocity profiles significantly de- 
part from the linear temporal behavior). A lower expansion rate 
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is a natural consequence of the compression induced by the over- 
taking flow. 

More surprising, some perturbed MCs are found to expand 
faster (larger £) than non-perturbed MCs. We conclude that such 
MC are probably observed after the main interaction phase with 
the overtaking flow, so that they expand faster than usual in order 
to sustain an approximate pressure balance with the surrounding 
SW (Fig. [3}. More precisely, as the overtaking flow disappears 
from the back of the MC, the MC is expected to tend towards the 
size that it would have reached without the interaction with the 
fast stream. Since it was compressed, it is expanding faster than 
usual to return to its expected size in a normal SW. 

5. Conclusions 

Our present results confirm and extend our previous work on 
the expansion of MCs. The non-dimensional expansion factor 
£ (Eq. [8} gives a precise measure of the expanding state of a 
MC. In particular, it removes the size effect which could give the 
false impression that parcels of fluid in large MCs are expanding 
faster. The value of f in non-significantly perturbed MCs is clus- 
tered in a narrow range, independent of their magnetic structure, 
but determined by the pressure gradient of the surrounding S W. 
The mean value of £ is also nearly the exponent of the solar dis- 
tance for the MC size (determined from the analysis of MCs at 
different solar distances, section 12.5) . 

However, for perturbed MCs, f has a much broader range, a 
result linked to its proportionality to the local expansion veloc- 
ity. So for perturbed MCs, ( is a measure of the local expansion 
rate and of the importance of the overtaking stream (i.e., a quan- 
tification of the influence of the MC/stream interaction on the 
expansion of the MC). 

Finally for non-perturbed MCs, £ is independent of the solar 
distance in the inner heliosphere. Presently we do not know how 
far this result extends to larger distances, even if it is an expected 
result as long as the flux ropes still exist. This will be the subject 
of a future study. 
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Fig. 5. Perturbed and non-perturbed MCs have a remarkably dif- 
ferent behavior of £ when they are plotted as a function of AV X . 
The drawing convention is the same as in Figs. I3I4I 
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Fig. 6. Cartoon of a possible evolution of the size of the MC with 
the helio-distance, showing the expected global expansion (thick 
solid line), an example of a non-perturbed MC (thin solid) and a 
perturbed MC (dashed line). 
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Fig. 7. Log-log plots of the SW total pressure and its compo- 
nents as a function of the solar distance. The solid line shows 
the least square straight line fitted to the total pressure as com- 
puted from the points marked with circles. The total pressure in 
the SW decreases as P SW (D) = P D~ 29 . 
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Fig. 1. Examples of two analyzed MCs that are not significantly perturbed by a fast flow. The MC center was observed at 07-Jan- 
1975 10:39 and 04-Mar-1975 21:37 UT, for panels (a,c) and (b,d), respectively. The vertical dashed lines define the MC boundaries. 
(a,b) V x is the observed velocity component in the radial direction from the Sun, expressed in km per second. The straight line is 
the linear least square fit of the velocity in the time interval where an almost linear trend is present (where the observations are 
presented as a solid line). The linear fitting is extrapolated to the borders of the MC, which are marked with circles. (c,d) B y is the 
magnetic field component, in nT, both orthogonal to the MC axis and to the spacecraft trajectory, while the solid line represents F y , 
which is the accumulated flux of B y (Eq. |7}. The extremum of F y (proxy of the cloud center) is indicated with diamonds. 
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Fig. 2. Examples of two analyzed MCs that are perturbed by a fast flow as seen in the upper panels. The MC center was observed 
at 30-Jan-1977 03:18 and 23-Jun-1980 12:25 UT, for panels (a,c) and (b,d), respectively. The external perturbations enter in a 
significant part of the MCs. The same quantities, as in Fig. Q] are shown. 
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Fig. 3. The panels a-d show the correlation analysis between proxies for MC expansion with different physical quantities. MCs are 
separated in two groups: perturbed (empty square symbol) and non-perturbed (filled circle symbol). The straight lines are the result 
of a least square fit for perturbed (dashed line), non-perturbed (thin continuous line), and for both set of MCs in the list of events 
shown in Table Q] (thick continuous line). AV X is defined by Eq. ©, Af = f out - t m , V c is the velocity of the MC center (or at the 
closest approach distance), and D is the distance to the Sun. The fitted values and the obtained correlation coefficient (r) are included 
as insets, considering the different groups. More significant correlation is present for the non-perturbed cases. 
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Fig. 4. Panels a-d show the correlation analysis that tests for the dependence of the non-dimensional expansion factor £ (Eq.O as a 
function of other parameters. Perturbed and non-perturbed MCs are represented as in Fig. [3] S corrected to l au and < B > C otreeted to l au 
are normalized to 1 AU using the size and field strength dependence on the distance, according to the relationship given in Eq. (O. 



